ON MOCKENHOUPT'S CONJECTURE IN THE HARDY-LITTLEWOOD 

MAJORANT PROBLEM 

SANDOR KRENEDITS 

Abstract. The Hardy-Littlewood majorant problem has a positive answer only for expo- 
nents p which are even integers, while there are counterexamples for all p ^ 2N. Montgomery 
conjectured that even among the idempotent polynomials there must exist some counterex- 
'f^' ' amples, i.e. there exist some finite set of characters and some ± signs with which the signed 

character sum has larger p^^ norm than the idempotent obtained with all the signs chosen 

^^ ' + in the character sum. That conjecture was proved recently by Mockenhaupt and Schlag. 

^N , However, Mockenhaupt conjectured that even the classical 1 + e^'^*^ ± e^'^*('''+^)^ three- 

J-j ' term character sums, used for p = 3 and k — 1 already by Hardy and Littlewood, should 

work in this respect. That remained unproved, as the construction of Mockenhaupt and 
Schlag works with four-term idempotents. In our previous work we proved this conjecture 
for fc = 0, 1, 2, i.e. in the range < p < 6, p ^ 2N. 

Continuing this work here we demonstrate that even the fc = 3, 4 cases hold true. Several 
refinement in the technical features of our approach include improved fourth order quadra- 

^ ' ture formulae, finite estimation of C^/G (with G being the absolute value square function 

of an idempotent), valid even at a zero of G, and detailed error estimates of approximations 
of various derivatives in subintervals, chosen to have accelerated convergence due to smaller 

I-S5 , radius of the Taylor approximation. 
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^ ! 1. Introduction 

We denote, as usual, T := M/Z the one dimensional torus or circle group. Following 

K> I Hardy and Litlewood [7], / is said to be a majorant to g if |^| < /. Obviously, then / is 

H ' necessarily a positive definite function. The (upper) majorization property (with constant 

- ■ ■' 1) is the statement that whenever / e L^(T) is a majorant of (? G Lp{T), then \\g\\p < \\f\\p- 



Hardy and Littlewood proved this for all p E 2N - this being an easy consequence of the 
Parseval identity. On the other hand already Hardy and Littlewood observed that this fails 
for p = 3. Indeed, they took f = 1 + ei + e^ and g = 1 — ei + e^ (where here and the sequel 
we denote efc(x) := e{kx) and e(t) := e^'"**, as usual) and calculated that ||/||3 < WgWs- 

The failure of the majorization property for p ^ 2N was shown by Boas [3]. Boas' con- 
struction exploits complex Taylor series expansion around zero: for 2k < p < 2k + 2 the 
counterexample is provided by the polynomials f,g := 1 -|- rei ± r'^"^^efc+2, with r sufficiently 
small to make the effect of the first terms dominant over later, larger powers of r. 
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Utilizing Riesz products - an idea suggested to him by Y. Katznelson - Bachelis proved f2] 
the failure of the majorization property for any p ^ 2N even with arbitrarily large constants. 
That is, not even \\g\\p < Cp\\f\\p holds with some fixed constant C = Cp. 

Montgomery conjectured that the majorant property for p ^ 2N fails also if we restrict 
to idempotent majorants, see pT| p. 144]. (A suitable integrable function is idempotent if 
its convolution square is itself: that is, if its Fourier coefficients are either or 1.) This has 
been recently proved by Mockenhaupt and Schlag in [TO]. 

Theorem 1 (Mockenhaupt & Schlag). Let p > 2 and p ^ 2N, and let k > p/2 be 

arbitrary. Then for the trigonometric polynomials g := [1 + efc)(l — e^+i) and /:=(! + 
efc)(l + Cfc+i) we have \\g\\p > 



Oddly enough, the quite nice, constructive example is given with a four-term idempotent 
polynomial, although trinomials may seem simpler objects to study. Indeed, there is a 
considerable knowledge, even if usually for the maximum norm, on the space of trinomials, 
see e.g. [SI [121 [IS]- Note that three-term examples are the simplest we can ask for, as two- 
term polynomials can never exhibit failure of the majorization property. In the construction 
of Mockenhaupt and Schlag, however, the key role is played by the fact that the given 4-term 
idempotent is the product of two two-term idempotents, the p^^ power integral of which then 
can be expressed by the usual trigonometric and hyperbolic functions. So even if four terms 
is a bit more complicated, but the product form gives way to a manageable calculation. 

Nevertheless, one may feel that Boas' idea, i.e. the idea of cancelation in the {k + ly^ 
Fourier coefficients works even if r is not that small - perhaps even if r = 1. The difficulty 
here is that the binomial series expansion diverges, and we have no explicit way to control 
the interplay of the various terms occurring with the ± signed versions of our polynomials. 
But at least there is one instance, the case of p = 3, when all this is explicitly known: 
already Hardy and Littlewood [7] observed that failure of the majorant property for p = 3 
is exhibited already by the pair of idempotents 1 -|- ei ± 63. In fact, this idempotent example 
led Montgomery to express (in a vague form, however, see [11], p. 144) his conjecture on 
existence of idempotent counterexamples. 

There has been a number of attempts on the Montgomery problem. In particular, led 
by the examples of Hardy-Littlewood and Boas, Mockenhaupt [H] expressed his view that 
1 + ci ± ek+2, where 2k < p < 2k + 2, should provide a counterexample in the Hardy- 
Littlewood majorant problem, (at least for k = 1,2). So we are to discuss the following 
reasonably documented conjecture. 

Conjecture 2. Let 2k <p <2k + 2, where fc G N arbitrary. Then the three-term idempotent 
polynomial Pk := 1 + ei + ek+2 has smaller p-norm than Qk '-=1 + ei — ek+2- 

We have proved this for fc = 0, 1, 2 in [8j. 

One motivation for us was the recent paper of Bonami and Revesz [4j, who used suitable 
idempotent polynomials as the base of their construction, via Riesz kernels, of highly con- 
centrated ones in L'^iT) for any p > 0. These key idempotents of Bonami and Revesz had 
special properties, related closely to the Hardy-Littlewood majorant problem. For details 
we refer to [4j . For the history and relevance of this closely related problem of idempotent 
polynomial concentration in L^ see [H |5], the detailed introduction of [S], the survey paper 
[T], and the references therein. 

As already hinted by Mockenhaupt's thesis [9], proving that 1 + ei ±6^+2 would be a coun- 
terexamle in the Hardy-Littlewood majorant problem may require some numerical analysis 
as well. However, we designed a way to accomplish this differently than suggested by Mock- 
enhaupt, for we don't know how to get it done along the lines hinted by him. Instead, in 
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[H] we used function calculus and support our analysis by numerical integration and error 
estimates where necessary. 

These methods are getting computationally more and more involved when k is getting 
larger. Striving for a worst-case error bound in the usual Riemann numerical integration 
formula forces us to consider larger and larger step numbers (smaller and smaller step sizes) 
in the division of the interval [0, 1/2], where a numerical integration is to be executed. There- 
fore, for k getting larger, we can as well expect the step numbers increase to a numerically 
extraneous amount, where calculations loose liability in view of the possibly accumulating 
small errors of the computation of the operations and regular function values - powers, log- 
arithms and trigonometrical or exponential functions - involved. Any reader would readily 
accept a proof, which with a certain precise error estimate refers to a numerical integration 
formula on say a few hundred nodes, but perhaps no reader would be fully convinced reading 
that a numerical tabulation and integration on several tens of thousands of function values 
led to the numerical result. Correspondingly, in this paper we settle with the goal of keep- 
ing any numerical integration, i.e quadrature, under the step number (or number of nodes, 
division number) N = 500, that is step size h = 0.001. 

Calculation of trigonometrical and exponential functions, as well as powers and logarithms, 
when within the numerical stability range of these functions (that is, when the variables of 
taking negative powers or logarithms is well separated from zero) are done by mathematical 
function subroutines of usual Microsoft excel spreadsheet, which computes the mathematical 
functions with 15 significant digits of precision. Although we do not detail the estimates of 
the computational error of applying spreadsheets and functions from Microsoft Excel tables, 
it is clear that under this step number size our calculations are reliable well within the error 
bounds. For a more detailed error analysis of that sort, which similarly applies here, too, 
see our previous work [S], in particular footnote 3 on page 141 and the discussion around 
formula (22). 

In view of the above considerations, instead of pushing forward exactly the same numerical 
analysis as done in [S] for k = 1,2 also for higher values of k, (which could have been done at 
least for some k, though), here we renew the approach and invoke a number of new features 
of the numerical analysis. These "tricks" will enable us to keep A^ below 500, and thus keep 
the invoked numerical calculations of quadratures reliable. 

First, instead of the classical and simplest numerical integration by using "brute force" 
Riemann sums, we apply a more involved quadrature formula (TTTI) . derived from Taylor 
approximation, which in turn allows us to keep the step number under good control. Here 
instead of the most famous Simpson rule, which uses only function values, we prefer a some- 
what more involved quadrature, calculating the approximate value of the integral by means 
of using also the values of the second derivative of the integrand. The gain is considerable 
even if not in order, but in the constant of the error formula. 

Second, as already suggested in the conclusion of [J8j, we apply Taylor series expansion 
at more points than just at the midpoint to := ^ + 1/2 of the t-interval {k,k + 1). This 
reduces the size of powers oi {t — to) , from powers of 1/2 to powers of smaller radii. The 
Taylor polynomial of degree 7, considered in [8J, had error size 2~^ due to the contribution 
of 1^ — to I ^ in the Lagrange remainder term, while here for fc = 4 the division of the t-interval 
to (4,4.5) and (4.5,5) results in 0(4~") in the respective error contribution. 

2. Notations and a few general formula for the numerical analysis 

Let A; G N be fixed. (Actually we will work with k = 3 or k = A only.) To set the 
framework, here we briefly sketch the general scheme of our argument, and exhibit a number 
of general formulae for later use in the analysis. 
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In the sequel we write F±{x) := 1 + e{x) ± e{{k + 2)x) and consider the p^^ power integrals 
f±{p) '■= /o \F±{x)\Pdx as well as their difference 



A{p) := f^p) - U{p) :- 



\F(x) 



\F+(xWdx. 



Our goal is to prove Conjecture [21 that is A{p) > for all p G {2k, 2k + 2). 

Let us introduce a few further notations. We will write t := p/2 G [k,k + 1] and put 



(1) 



(2) 



G+(x) 



F+(x) 



1/2 

9±{t):=\h{2t) = j G 



iXj j HjJu • 



1/2 

d{t):=\/\{2t)=g^{t)-g^{t) = j [G{ 



\x 



G\{x)\ dx. 



So we are to prove that d{t) > for k < t < k + 1. First we derive that at the endpoints d 
vanishes; and, for later use, we also compute some higher order integrals of G±. 



Lemma 3. Let p E N with 1 < p < k + 1. Then we have 
(3) 



G^ 



l^±l' 



p-{k+2) 
iy=0 



with 
(4) 

where p :- 
(5) 



a±{iy) ■■= (±1) 






p- p 
A 



u 



k + 2 



and A := u — p{k + 2) is the reduced residue of u mod k + 2. Therefore, 

^/^ p-ik+2) 






/n particular, jl''^ \G+\p = J^^"^ \G-\p for all < p < k + 1 and thus d{k) = d{k + 1) = 0. 

Remark 4. By similar calculations one can compute a±{i') even for higher values of p as 
well. E.g. m the range k+2 < p < 2k+3 we have a±{u) = (±1)^ | (|;) (p-^) ± (^^ J {x+(l,Xl-)) ] ■ 
That we will not use, however. 

Proof. In the trinomial development of (l+Ci ±6^+2)^ the general term coming from choosing 
a times ±6^+2 and r times ei (and then necessarily p ~ a — t times the constant term 1) 
has the form (±1)'^(^) (^~'^)- This to contribute to a±(z/) we must have z/ = a{k + 2) + r, a 
condition which forces u = r mod k + 2. Now if p < fc + 1, we also have < r < fc + 1, 
so the number r of choosing eis is exactly A, the mod k + 2 reduced residue of u, and 
consequently a = (z/ — A)/(fc + 2) = p. That results in formula (jl]) for a±(i/), while ([5]) 
follows by Parseval's formula. Whence the assertion is proved. D 
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To start the analysis of G{x) := G±{x), let us compute its x-derivatives. We find 
G±{x) = 3 + 2{cos(27r2;) ± cos((2A; + 2)ttx) ± cos((2A; + A)7ix)}, 

(6) G'^{x) = -Att sin(27ra;) T (4A; + 4)7r sin((2A; + 2)7ra;) T (4A; + 8)7r sin((2A; + 4)7rx)), 
GI{x) = -8n^ cos(27rx) T 8(A; + 1) V^ cos((2A; + 2)nx) T 8(A; + 2)V cos((2A; + 4)7rx)}, 

and in general 

■ {sin(27rx) ± (A; + 1)^™+^ sin((2A; + 2)7rx) ± (A; + 2)^"^+^ sin((2A; + 4)7rx) } , 

G^l"'\x) = 2{-Arn^"' 

■ {cos(27ra;) ± {k + 1)^" cos((2fc + 2)nx) ± {k + 2)^"^ cos((2fc + 4)7ra;)} . 

Consequently we have 

(7) ||G'±|U<9=:Mo, 

||Gi"^^||oo < 2'"+V™{l + {k + ir + {k + 2)™} =: M^k) =: M„ (m = 1, 2, . . . ). 

We encounter a new phenomenon, compared to [8j|, when k = 3, since here G'+(a;) does 
not have a positive lower bound: we in fact have G+{l/3) = 0. (Let us note in passing that 
for G_ we have miuT G_ ^ 0.282... > 1/4 - but we do not use this in the following.) 

For higher x-derivatives of the composite functions G"^log''G'+, needed in our analysis, 
vanishing of G^ causes concerns for occurring negative powers of (?+ after differentiation, 
while the appearance of log-' G+ invoke concerns of blowing up calculations and estimates in 
view of "logO = oo". The first problem we resolve by a comparison of G+ to G'^, always 
present in the numerator, while the second difficulty will be taken care of by using only 
continuous functions u"log''w, with a > 0, 6 > 0, of f = G+{x). Although all this can be 
avoided, when G_ is strictly bounded away from zero, for a possibly better estimation we still 
calculate the same comparative estimates even for G-. (Similarly, the idea of comparison of 
G'^ and G± could be used for higher k as well, whether or not the functions G± vanish.) 

So we want to compare G' and vG = IF I, more precisely G'^ and G. Note that G' = 

m ■ m' ^ 2VG- {VaJ. Anothe. heu Jicd .ea=o.ng toiustify the .a.ch fo. a bound of 

G'^/G, is that G > 0, hence whenever G = we necessarily have G' = 0, and the multiplicity 
m of any zero of G being an integer (as G is an entire function), we conclude m >2: so G'^ 
has a zero of order 2{m — 1) > m. 

So we start the search for a bound on G'^/G±. To this end we write u = cosf with 
V = 271X and calculate 

G'iix) = (47r)2(sinw ± (A; + 1) sin(A: + l)v ± {k + 2) sm{k + 2)vY 

(8) = 167r'(l - u^) [l±{k + l)Uk{u) ±{k + 2)Uk+i{uf , 

1 TT / \ sin((m + l)v) \ . , i ^, i i i • i r i 

where UmVu) := (v := arccosMj is the m-th Chebyshev polynomial of the 

sint> 
second kind. 

We are to compare this and 

G^(x) = 3 + 2 cost; ± 2 cos(A; + l)t; ± 2 cos(A; + 2)v = 3 + 2u± 2Tk+i{u) ± 2Tk+2{u), 

where here Tm{u) = cos{mv) {v := arccos-u) is the mth Chebyshev polynomial of the first 
kind. 
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In all, G'^ /G is always an entire function of x, and substituting u 
have the formula 



6 
cosw = cos27rx we 



(9) 



G'^{x) _ 167r^(l - u^) [l±{k + l)Uk{u) ± (A: + 2)Uk+i{u)\ 
G{x) 3 + 2M±2Tfc+i(n)±2Tfc+2(n) 



In the paper [8] we used Riemann sums and the standard Riemann sums approximation 



formula 



^1/2 



tions of the form $ 



1 V^^$ 



2N 



ra-1/2 
2iV 



< 



II*"!!. 



when numerically integrating func- 



192Ar2 ' 

H := G* log-' G along the x values. 

A new feature of the present approach is that for better approximation we now improve 
the numerical integration method by means of invoking a quadrature formula. This was 
not feasible for small t, as higher derivatives of the composite function H lead to G in the 
denominator: the im}^ derivative in general results in the occurrence of G*~™, and negative 
powers of G bear the risk of blowing up all of our estimates. This can be remedied a little 
by comparison of G'^ to G, a lucky possibility explained above. This was already utilized 
in |H] to control 2"^^^ derivatives of H, and we'll make use of it here, too for /c = 3, when for 
some integrals (some occurring H functions) t can be as small as t = 3, while we need to 
control 4*^ derivatives of H in view of error terms of the quadrature formula we use. With 
this additional consideration the 4*^ derivatives of H can always be controlled for A; > 3. 
(For /c = 0, 1, 2, settled in |i8j, this could not have been possible.) 

For remaining self-contained, we deduce here the otherwise well-known quadrature formula 
what we want to apply. This starts with the 3"^^^ order Taylor polynomial approximation (with 
the so-called Lagrange error term), valid for four times continuously differentiable functions ip: 



ip{x) = ip{xo) + ip'{xo)ix - Xo) + 



2 ' ^' 6 

Integrating over a symmetric interval [xq — q,xo + q] leads to 



^"(^°) (X - xof + ^^(x - xof + ^'""^^--^ - ^ ^^ 



24 



la; - Xo) 



xo+q 



Q 

ip{x)dx - { ip{xo)2q + ip"{xo) — 



po-q 



XQ+q 



^'""iC. 



24 



(x — xoVdx 



po-g 

xo+q 



JV I 



< max \ip [x) 

[xo-q,xo+q] 



(x-Xq)^ jy q^ 
ax < max \(p [x)\ — . 

24 [xo-q,xo+q] 60 



xo-q 



Applying the same formula for A^ intervals of the form [x„ — /i/2, x„ + h/2], where h 
{b — a)/N and x„ = (n — l/2)h + a with n = 1, . . . , A^, we obtain 



(10) 



/ ^ ~ 5Z 1 V{Xn)h + v"{Xn)— \ 
{ n=l ^ ^ 



h' 



N 



max \(I> (x) < r 09 loo- 

-- „=il— nl<t " fin 95 11^ II- 



eo^s 



This leads to the following quadrature formula^ 



Note the noticeably better error estimate, not in order but in constant, than one would obtain by more 
customary Simpson type rules. This is due to the use of second derivatives, which in our case will still be 
calculable, explicit formulae. 
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Lemma 5. Let (f be a four times continuously differentiable function on [0,1/2]. Then we 
have 



:iii 



1/2 N 

/ ip{x)dx - Y^ 

•I n=l 



V 



2n 



AN 



1 

2iV 



V 



2n-l 



1 



AN / 192A^3 



< 



60 2^m^' 



Let us start analyzing the functions 
(12) H{x):= Htj,±{x) := G*±(x) log^' G±{x) 



[xG [0,1/2]) (te[k,k + l],jEN). 



To find the maximum norm of Htj^±, we in fact look for the maximum of an expression of 
the form f*| logf p, where v = G{x) ranges from zero (or, if G 7^ 0, from some positive lower 
bound) up to UGHoo < 9. For that, a direct calculus provides the following. 

Lemma 6. For any s > and m G N the function a{v) := as^miv) '■= f*| logf I™ behaves on 
[0, 00) the following way. It is nonnegative, continuous, continuously differentiable, (apart 
from possibly in case s < 1), has precisely two zeroes at and 1, and it has one single 
critical point vq = exp(— m/s). Consequently, it has exactly one local maximum point at Vq 
where its local maximum is (— ) , furthermore, the function increases in [0,fo] and also on 
[1,00), and decreases on [vq,!]. Therefore for any finite interval [a,b] C [0, 00) we have 



(13) 



raaxaiv) 

[a,b] 



if a < b < Vq, 
if a < Vq < b < 1, 
ifa<vo, 1 < b, 
if Vq < a < b < 1, 
if Vq < a < 1 < b, 
ifl<a<b. 



a{b) 
a{vQ) 

max{Q;(fo), «(&)} 
a{a) 

max{a(a), a{b)} 
(aib) 

In particular for [a,b] = [0,9] we always have 

'i^r ifmls>^ 

\es J J I (jQ 

9Mog™9 ifm/s<^, 
where (Tq ~ 0.126... is the unique root of the equation cr9°" = l/(elog9). 

For the application of the above quadrature (ITT|) we calculate (c.f. also [8J) 

H"{x) := H'/.^^{x) = G"{x)G''\x) log^-^G{x) {tlogG(x) + j} 

(15) + G'\x)G'~\x) log^-2 G{x) {t{t - 1) log^ G{x) + j{2t - 1) logG(x) + j{j - 1)} . 

However, the error estimation in the above explained quadrature approach forces us to 
consider even fourth x-derivatives oi H = Htj^±. In order to calculate 

4 



(14) 



a. 



max aiv) 

[0,9] 



max <; — 

\es 



9Mog'"9} 



(16) H 

we start with computing 



IV 



Tn=0 



m 



{G'y"'>{\og'G) 



(4— m) 



(17) 



(G*)' = tG'-'G' 

(G*)" = t{t-l)G'-''G'^ + tG'-^G" 

(G*)'" = t(t-l)(t-2)G*-3G'3 + 3t(t-l)G*-2G'G" + tG*-iG'" 

(Gy = t{t - l){t - 2){t - 3)G'~^G'^ + 6t{t ~ l){t - 2)G'-''G'^G" 

+ 3t{t - 1)G'-^G"^ + At{t - l)G'-^G'G"' + tG'-'G^^ 
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and denoting L := logG also 

{Ur = —J [U - l)(j - 2)L^-' - 3(j - IW-' + 2U-'] 
(18) + -^3j[(j - \)V-^ - U-'] + -^3L^-\ 

(L^y^ = ^3 [[{j - l)(j - 2)(j - 3)L^-^ - 6(j - l)(j - 2)U-' + ll(j - l)^-' - QU'^] 
+ ^6j [(j - l)(j - 2)U-' - 3(j - l)U-' + 2L^-i] 

Inserting flT7|) and flTSj) into (1161) we arrive at the desired general formula for H^J^ as follows 

H^y = G'-'G''{jij - l)(j - 2)(j - 3)L^~^ + [At - 6]j(j - l)(j - 2)L^~' 

+ [et^ - 18t + ll]j(j - l)L^'"2 + [2t^ - 9t^ + lit - 3]2jL^-i + t(t - l)(t - 2)(t - 3)lA 

+ 6 • G'-'G''G"{j{j - l)(j - 2)U-' + 3(t - l)j(j - 1)^^^-' + [3t^ - 6t + 2)]jL^-i 
(19) 

+ t(t - l)(t - 2)l4 + 4 ■ G*-'G"G""|j(j - 1)^'"' + (2t - l)jX^-i + t{t - 1)l4 

+ G*-^G'^ {jL^~' +tL^}+3- G'~^G"^lj{j - l)U-^ + (2t - l)jX^-i + t{t - 1)l4. 

At all occurrences we will need an estimate for ||if ^^||oo in order to apply it in the numerical 
quadrature formula. Therefore, we now start estimating the above expression. For a shorter 
notation we write v := G{x) G [0,9] and d, := \L\ = \ logf |. As a first step we thus find for 
j = 1, 2, 3, ..., t > 3 the estimates 

in'^^ix)] < v'-'Mt{j{j - l)(j - 2)(j - 3)f-' + [At - 6]j(j - l)(j - 2)f-' 

+ [et^ - 18t + ll]j(j - 1)^^'"^ + [2t^ - Qf^ + lit - 3]2j&-'^ + t{t - l)(t - 2)(t - 3)^4 

+ 6 ■ v'-^MlM2[3{j - l)(j - 2)P-^ + 3(t - l)j(j - l)P-^ + [3t2 - 6t + 2\je-^ 
(20) 

+ t{t - l)(t - 2)P\ + t;*~^M4 {jF~^ + t£^} 

+ t;*-2(3M2 + 4MiM3)|j(j - 1)^^'"^ + (2t - l)]^-^ + t(t - l)e\. 



Furthermore, to be used typically for smaller values oi v = G{x), that is to say only for 
< w < 3, we can derive a different estimation whenever some constant M* := M*{k) is 
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known satisfying HG'^/GHoo < M*. Namely, we then have 

\H'''{x)\ < v'-^M*^[j{j - l)(j - 2)(j - ?,)&-' + [At - 6]j(j - l)(j - 2^-' 

+ [et^ - 18t + ll]j(j - 1)^^"^ + [2t^ - 9t^ + Hi - 3]2j£^-^ + t(t - l)(t - 2)(t - Z)P\ 
+ 6 ■ t;*-2M*M2{j(j - l)(j - 2)F"^ + 3(t - l)j(j - l)P'^ + [St^ - 6t + 2)]jF-i 



+ t(t - l)(t - 2)F + 4 • v'-^-^VWmAj{j - l)f-^ + (2t - 1)3 f-^ + t{t - 1)F 



f211 



+ v*-iM4 {j^^-^ + W]+?,- w*-2m||j(j - 1)^^"^ + (2t - l)jF^^ + t(t - 1)£A. 



Furthermore, estimating by means of A := max(£, 1) and using l^ ^,P ^, P ^, P '^ < A-' and 
also 2^/v < 1 + 1" we are led to 

\H'''{x)\ < v'-'A^ (m*'[j{j - l)(j - 2)(j - 3) + [At - 6]j(j - l)(j - 2) 

+ [6t^ - 18t + ll]j(j - 1) + [2t^ - 9t^ + Hi - 3]2j + t{t - l)(t - 2)(t - 3)} 
+ 6 ■ M*M2[j{j - l)(j - 2) + 3(t - l)j{j - 1) + [3t2 - 6t + 2)]j + t{t - l)(t - 2)} 

+ 2VIFMs{jij - 1) + (2t - l)j + t(t - 1)} + 3M^{jij - 1) + (2t - l)j + tit - 1)} 
(22) 

+ v'-'A^h- v/M^MsJjXj - 1) + (2t - l)j + t(t - 1)} + M4 {j + t} J . 

3. The proof of the k = 3 case of Conjecture [2] 

When fc = 3, let us start with a few concrete numerical estimates of the functions G^ 
and H^^\ For fc = 3 we need U-^iu) = Au{2u'^ - 1), U^iu) = 16m^ - 12^2 + 1 and n{u) = 
8u^ -8^2 + 1, T5{u) = 16u^ - 20u^ + 5m. Writing these in (jO]) yields for G+ 

G'l , , 167r2(l - u^) [80m^ + 32^3 - eOw^ - I6n + 61^ 



G^ ' ' 32^5 + 16^4 _ 40u3 - 16^2 + 12m + 5 

_ 167r2(l - m2)[40m3 - 4^2 _ 28m + 6]^ 
^ ' ~ 8n3 - 4^2 - 8n + 5 ' 

canceling the common factors of (2-u + 1)2. Note that the denominator is now non-vanishing 



in the interval [—1,1], as its minimum is ~ 0.12, attained at ~ 0.76759.... Thus 

the above rational function can be maximized numerically on the range m G [—1,1] of m = 
cos(27rx), the maximum being f« 3699, so 

(24) G'l{x) < 3700G+(x). 

Although G- does not vanish, for a possibly better estimation for small values of G^{x), 
we still work out a bound on G'^/G-. Again with u = cosf and v = 2ttx we get from (j9]) 

(25) ma. ^(x) = ,„a. l^^d - .') |80.- + 32.3 - 60„. - 16. + 4f ^ 

^ ' [-1,1] G- ^ ' [-1,1] -32w' - 16^4 + 40m3 + 16^2 - 8n + 1 
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Note that now the denominator does not vanish and there is no singularity to make the 
numerical maximization difficult. Summing up, we find 



(26) 



G'l 



< M*(3) := 3900. 



G± 

The next step is, as in [8j, to see that d^^\3) > for the first few values of j = 1, 2. 

Lemma 7. We have d'{3) > 0. 

Remark 8. A preliminary numerical calculation yields the approximate value d'{3) ~ 0.01401... 
We don't need the concrete value, but this information suggests us the proper choice of the 
targeted error bound of 6 = 0.007 below. 

Proof. From ([2]) we clearly have 

1/2 1/2 

rfO)(t) = g^^\t) - g^l\t) = I G'_{x) log^' G^{x)dx - j G\{x) log^' G+{x)dx 



1/2 1/2 

(27) = / Htj-{x)dx — / Htj^+{x)dx. 



Now we calculate the value - that is, these two integrals - numerically for t = A; = 3 and 
j = 1. Both integrals should be computed within the error bound 5 := 0.007. Invoking 
Lemma [5] we are left with the estimation of H-f^sYitlloo- The general formula of flT9|) now 
specializes to 

(28) 

H^^ = 6^ + G'^G'\m + 36L) + GG"G'"(20 + 24L) + GG"^{lh + 18L) + G^G^^il + 3L). 
G 

We now estimate \H^^ {x)\ distinguishing two cases, the first being when v := G{x) > 3. 
Inserting the estimates of HC^™) ||oo from (jTj) for m = 0, 1, 2, 3, 4, we get from 



\H^^{x)\ < 6^^^ + (407r)2(8 • n^ ■ 42) (66 + 36 log t;) + t;(407r)(167r3 ■ 186) (20 + 24 log w) 

V 

+ f (871^42)2(15 + 18 log t;) + t;2(327r^882)(l + 3 log t;) 

(29) = 7r^{15, 360, 000/t; + 35, 481, 600 + 19, 353, 600 log w + 4, 074, 240t; 

+ 4, 889, 088vlog v + 28, 224^;^ + 84, 6721;^ log v}, 

which is clearly an increasing function of f = G{x) for v > 2, e.g. Therefore substituting the 
maximal possible value f = 9 we obtain in this case 

(30) \H^^ix) I < 22, 444, 818, 695 < 2.3 ■ 10^°. 

For smaller values of f = G{x) we estimate (1281) the same way as it is done in general in 
( ETj) . with Mm in ([7]) and M* in (^ (or, we substitute t = fc = 3 and j = 1 in (12ID and use 
the numerical values of M^ and M* as said) . This yields 

\H^^{x)\ < 1.2 ■ 10% + (6.7t;logt; + 1.2^;^/^ + 1.4t;=^/2iQg^) . ^gS + (2.8t;2 + 8.3t;2logt;) ■ 10^ 

(Also, we could have substituted t = k = 3 and j = 1 in f l2^ and apply ([7]) and f EBjl in 
that.) The function on the right hand side takes its maximum on [0,3] at f = 3, thus 

\H^^{x)\ < 7.014 ■ 10^ < 8 ■ 10^ 
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is obtained in this second case. In all, we find ||if^^(x)|| < 2.3 • 10^°, hence in the numerical 
quadrature formula fill I) the error is estimated by 

2.3-101° 



60 • 210 Ar4 
To bring this down below 6 = 0.007, we need to chose the step number A^ as large as to have 



2-3 -101° , . .r nr J 2.3-101° 

<6 i.e. N>Nn:=\ ^ 8Q.... 



60 ■ 210 A^4 - V 60 ■ 210 ■ 0.007 

Calculating the quadrature formula with A^ = 100, we obtain the approximate value 0.014012641..., 
whence d'{3) > 0.014012641... - 2 ■ 0.007 > 0. D 

Lemma 9. We have d"{3) > 0. 

Remark 10. Preliminary numerical calculation yields d"{3) ~ 0.087602... . 

Proof. From (1271) now we calculate the value - that is, these two integrals - numerically for 
t = k = 3 and j = 2. Both integrals should be computed within the error bound 6 := 0.04. 
As before, invoking Lemma [5] we are left with the estimation of ||-f^3 2^±||oo- The general 
formula of (IT9|) now specializes to 

H^^ = G'^G"{36L'^ + 132L + 72) + GG"\l8L^ + 30L + 6) 
(31) + GG'G"'{2AL^ + AOL + 8) + G^G^^{3L^ + 2L) + ^(12L + 22). 

Gr 

Similarly as before, in the estimation of \H^'^{x)\ we distinguish two cases. Namely, we 
separate cases according to v := G{x) > e or < f < e. For the case when e < f < 9, i.e. 
1 < L < log 9, application of ([7]) after substituting t = k = 3 and j = 2 in (l20l) (in other 
words, using ([7]) in estimating ([SI])) yields 



\H^^{x)\ < (407r)23367r2(36L2 + 132L + 72) + G{3367cy{18L^ + 30L + 6) + G407r29767r3 

tOTT 

'g' 



(AOtt)'^ 
■ (24^2 + 40L + 8) + ^^28, 22471^(3^^ + 2L) + ^ ^ (12L + 22) 



1 OQ 81 fi T^Q 

(32) < 7r^{288, 064, 517.... + 43, 137, 255. ..t; + 344, 030.. .t;^ + }, 

V 

which is, by easy calculus, an increasing function oi v = G{x) for v > 2, e.g. Therefore 
substituting the maximal possible value v = 9 yields in this case 

(33) |iJ^^(x)| <7-10i°. 

For smaller values of G{x) when < i; := G{x) < e, in fl22|) we substitute t = k = 3 and 
j = 2 and then use ([7]) for m > 2 and also fl2B]) . leading to 

\H^^{x)\ < vA'^ (34:M*^ + 240M*M2 + 3QVm*M3 + 5AM^) + v^A^ (36^/M*M3 + 5mA 

(34) < e (34 • 3900^ + 240 ■ 3900 ■ 3367r2 + 36^3900 ■ 304071^ + 54 ■ 336^71^) 

+ e^ (3673900 ■ 30407r3 + 5-28, 2247r^) ^ 13, 700, 830, 408 < 2 ■ 10^ 

applying also that on [0, e] f A^ < e and t>^A^ < e^. 
Summing up, 

II ff^^ll ^7 inio 



1IO 
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hence in the numerical quadrature formula (11 II) the error is estimated by 

7 • 10^0 



60-2iOAr4- 

We need to chose the step number A^ large enough to bring this error below 6 = 0.04, i.e. to 
have 

I 7 . inio 

N> No:= \ -7^ 7 ^ 73.05.... 

V 60 ■ 210 ■ 0.04 

Calculating the quadrature formula with A^ = 100, i.e. step size h = 0.005, we obtain the 
numerical approximate value 0.08760174..., so d"{3) > 0.08760174... - 2 ■ 0.04 > 0. D 

Our next aim will be to show that d" is concave in [3,4], i.e. that d^^ < 0. That will 
be the content of Lemma fT2l To arrive at it, our approach will be a computation of some 
approximating polynomial, which is, apart from a possible slight and well controlled error, 
a Taylor polynomial of d'^ . 

Numerical tabulation of values gives that d'^ is decreasing from d^^{3) ~ —0.068447... 
to even more negative values as t increases from 3 to 4. Thus our goal is to set n G N and 
6j > 0, {j = 0, . . . ,n + 1) suitably so that in the Taylor expansion 

(35) d^''{t)=J2^-r^{t-l) +Rn{d^\t), Rn{d^''^t) = ^^^[t-lj 

the standard error estimate 



\Rn{d'^,t)\< 

(36) < 



< 



||-^g,n+5,+ ||Li[0,l/2] + \\H^,n+5,^\\L^[0,l/2] ^-(ri+l) 

{n + l)\ 

1 II W II _L 1 II W II 

2 ll-"C.™+5,+ lloo "T 2\\''^^,n+5,~\\oo 

(n + l)!2"+i 
max3<^<4 ||if^_„_|_5^+||oo + max3<^<4 ||if^_„_|_5^_||oo 



(n + l)!2"+2 

provides the appropriately small error \\Rn{d^^ , ■)||oo < ^n+i, while with appropriate approx- 
imation 5^- of rf(J+^)(7/2), 

\dU+^)(i)-dA 



(37) 



S^^'\l)~d, f^_7y 



-^, '-^<S, (j = 0,l,...,n). 



Naturally, we wish to chose n and the partial errors 6j so that X^^o ^j < S := 0.068, say, so 
that d^^{t) < Pn{t) + 5 with 

(38) ^«W-Ef(*-i)'- 

Here the approximate values dj will be obtained by numerical integration, using the quad- 
rature formula flTTl) in approximating d^^'^il /2), which has the integral representation fl27j) 
with j = 0, . . . , ra. To be precise, we apply the error formula of fITT]) with Nj G N steps, 
where Nj are set in function of a prescribed error of approximation r/^, which in turn will be 
set in function of the choice of 5j . 

So now we carry out this programme. First, as G±{x) G [0,9], \G%.{x)\o^ G±{x)\ < 
max[o,9] \v^ log™ v\ = a| „j, which we consider with ^ G [3, 4] and m = n + 5 > i. By Lemma 
[6] we derive for all n < 18 that 

(39) \\H^,n+5,±{x)\\oo < 9Mog"+^9 = 6561 ■ 2"+Mog"+^ 3 (3 < ^ < 4, 1 < n < 18). 
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In view of (ESD this yields \RJS^\t)\ < '- 1 < 0.011 =: ^n for n = 10. 

(n+ 1)! 

Now we must set 6q, . . . , 6io, too. So let now 6j = 0.005 for each j = 0, . . . , 10. The goal 
is that the termwise error f l37|) would not exceed 6j, which will be guaranteed by Nj step 
quadrature approximation of the two integrals defining S^^^\7/2) with prescribed error 77^ 
each. Therefore, we set rij := 5j2^j\/2, and note that in order to have (!37|) it suffices that 



(40) 



iv, > n; 



'-"7/2j+4,±ll«2 

60 ■ 21%- 



\H^^ II 

■ 2^0j\2^-^j 



according to Lemma [51 So at this point we estimate WH^Y^ •.^4 _|_||oo for j = 0, . . . , 10 to find 
appropriate values of N*. 

Lemma 11. For j = 0, . . . , 10 we have the numerical estimates of TableUlfor the values 
of \\HjY ._j_4_|_||oo- Setting 6j = 0.005 for j = 0, . . . , 10 the approximate quadrature of order 

500 =: A^ =: Nj > N* with the listed values of N* yield the approximate values dj as listed in 
TableUl admitting the error estimates (IHT]) forj = 0, . . . , 10. Furthermore, ||-Rio(c^^^)i^)||oo < 
0.011 =: ^11 and thus with the approximate Taylor polynomial Pio{t) defined in (1381) the 
approximation \d^^(t) — Pio{t)\ < 6 := 0.068 holds uniformly for t G [3,4]. 



Table 1. Estimates for values of 



with 5j := 0.005, and values of dj with A^ : = 



I H^^ II 

N,: 



corresponding values of A^* 
= 500 for j = 0, . . . , 10. 




1 

2 
3 

4 
5 
6 
7 
8 
9 
10 



estimate for 
3.3 
9.1 
2.5 
6.8 
1.9 
4.8 
2.8 
2.6 
2.7 
2.9 
3.4 



H^^ 1 

-'-'7/2j+4,±l °° 


A^ 


10^2 


383 


10^2 


415 


10^3 


378 


1013 


310 


10^4 


239 


10^4 


169 


10^5 


142 


lO^s 


128 


w 


115 


10^8 


101 


10^9 


88 



dj 
-8.097236891 
-37.59530251 
-141.3912224 
-468.2134571 
-1423.831595 
-4074.963995 
-11,148.7318 
-29,465.89339 
-75,792.43387 
-190,751.6522 
-471,634.7482 



Proof. We start with the numerical upper estimation of H, 



IV 

7/2,i,± 



[x) for 3 < a; < 4, where 



now in view of the shift of indices we need the estimation for 4 < j < 14. All what follows 
is not sensitive to j < 14, but it is convenient that j > 4, as otherwise in some derivatives 
the powers of L{x) = logG'(x) would diminish, changing the formula slightly. 

When V = G{x) > e, substitution of t = 7/2 in (12(1 while using the estimates ([7]) of 
IIG^™)]! and £ < log9 < 2.2 yields the estimate 

{x)\<2.2^{^ + B.Vg}, 



I fjiv 



where 



A:-- 
B :-- 



A, 



9 ■ 10^ [f - 2f + 105j2 + 134j + 154] 
2.8 • 10^j=^ + 3.9 ■ 10®j2 + 5.7 ■ lO^j + 1.2 



10 



10 
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This last function is a strictly convex function oi u := y/v = \jG{x) - so it must have 
a unique minimum and two monotonic parts before and after the minimum point. Easy 
calculus yields that the minimum is located at fo := -g. Now vq is less than e, when j < 14, 
hence then the function is increasing on [e, 9] and it achieves its maximum at f = 9. When 
j = 14, the minimum falls inside [e, 9], and the maximum is the maximum of the values at e 
and 9, but the latter being much larger, we again find that our estimate is maximized taking 
V = 9. Finally, substitution of f = 9 yields 

(41) \Hl^j2,j,±{x)\ < 2.V ■ {3 ■ lO^j^ + 7.8 ■ lO^j^ + 1.5 ■ lO^j^ + 1.8 ■ 10^°j + 3.7 • 10^°} . 

When G{x) < e, we substitute k = 3 and t = 7/2 in fl22l) . and apart from the values of 
M2, M3 and M4 we also use the last estimate of fl26p . Further, we write in v^^'^A^ < e-f^/^A-' 
and, by means of Lemma [6l max[o,e] f^/^A-' = max (e^^^, (3^)"') which then yields 

(42) 

l^7Ai,±(^)l < max [ e3/^ ( — J j {1.6 ■ lO'^f + 1.1 ■ 10^ f + 5.6 ■ 10^ f + 1.6 ■ lO^j + 1.9 ■ 10^} 

From here we take the maximum of (14T]) and the above (H2|) for all j = 6, . . . , 14, which 
means using ( HT]) up to j = 9 and then (H2]) for j = 10, . . . , 14, leading to the upper estimates 
of \\Hy^ ,±||oo as hsted in Tabled] 

Finally, we collect also the resulting numerical estimates of N* - as given by the formulae 
( HOj) - in Table [Hand furthermore list the accordingly computed values of dj, too, applying 
the numerical quadrature formula (ITT1) with step size h = 0.001, i.e. N = Nj = 500 steps. D 

Lemma 12. We have S^{t) < for all 3 < t < 4. 

Proof. We approximate d^^{t) by the polynomial Pio(i) constructed in fl55]) as the approx- 
imate value of the order 10 Taylor polynomial of d^^ around to •= 7/2. As the error is at 
most (5, it suffices to show that pit) := Pwif) + 5 < in [3,4]. 
Now Pio(3) = -0.068458667... so Pio(3) + 5 <Q. 

Moreover, p\t) = P{oW = E!1i / ■ \ m (^ - 7/2)^-i and p'{3) = -4.00969183 < 0. 

(j - 1)! 
From the explicit formula oi p{t) we consecutively compute also p"{3) = —23.12291565 < 0, 

p"'(3) = -93.80789264 < andp(^)(3) = -324.0046433, p^^\3) = -978.7532737..., p''^\3) = 
-3144.062078..., p(^)(3) = -5587.909055..., all < 0. Thus p'-J\3) < for j = 0, . . . , 7. 

Therefore in order to conclude p{t) < for 3 < t < 4 it suffices to show that p^^\t) = 
ds + dg{t — 7/2) + ((iio/2)(t — 7/2)^ stays negative in the interval [3, 4]. However, the leading 
coefficient of p'-^^ is negative, while it is easy to see that the discriminant A := dg — 2d^diQ of 
p'^^' is negative, too: A ?^ —3.511 ■ 10^°. Therefore, the whole parabola of the graph of p^^-* 
lies below the x-axis i.e. p^^\t) < (Vt G M). It follows that also pit) < for all t > 3. D 

Proof of the k = 3 case of Conjecture\^ Since d{3) = (i(4) = 0, and d'{3) > 0, d takes 
some positive values close to 3; so in view of Lagrange's (Rolle's) theorem, d' takes some 
negative values as well. Therefore, d' decreases from a positive value at 3 to some negative 
value somewhere later; it follows that d" takes some negative values in (3,4). Also, d" is 
concave and d"{3) > implies that d" changes from positive values towards negative ones; 
by concavity, there is a unique zero point r of d" in (3, 4), where d" has a definite sign change 
from positive to negative. 

It follows that d', starting with the positive value at 3, first increases, achieves a maximal 
positive value at r, and then it decreases, reaches zero and then eventually negative values. 
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as seen above. That is, when it becomes zero at some point cr, it already has a negative 
derivative, and it keeps decreasing from that point on. So d' is positive until a, when it has 
a strict sign change and becomes negative until 4. Therefore, d increases until a and then 
decreases till 4; so d forms a cap shape and it is minimal at the endpoints 3 and 4, where it 
vanishes. It follows that rf > in (3,4). 

This concludes the proof of the /c = 3 case of Conjecture [21 D 

4. The case A; = 4 of Conjecture [2] 

First of all let us record that in case /c = 4 in (j9]) we are to deal with G±{x) = 3 + 2m ± 
2\T^{u) + Tq{u)] {u = cos27ra;) so putting in T^i^u) = 16u^ — 20u^ + 5u and Tq{u) = 
32m® — 48u^ + 17u^ a numerical calculation of the occurring polynomials give 

(43) minG+ ^0.0946... and minG. ^ 0.02776... 

T T 

Therefore in case fc = 4 we can estimate i := \L\ = \ logG(x)| < | log(0.027)| < 3.7 for both 
signs of G±. 

Lemma 13. We have (i'(4) > 0. 

Remark 14. By numerical calculation, d'{4) ^ 0.0062067.... 
Proof. From dH]) with t = 4, j = 1 

(44) H^^ = G'\50 + 24L) + 6GG'^G"(26 + 24L) + 4G^G'G"'{7 + 12L) + G'^G^^{1 + AL) 

+ ?,G^G"^{1+12L). 

From this, d < 3.7 and the estimates ([7]) we get by plain substitution as before ||i7-'^^||oo < 
1.6-1012 

To bring the resulting error estimate down below 5 = 0.003, we need to chose the step 
number N as large as to have 



1.6-1012 / 1.6.1012 

<6 i.e. N> No:= \ L^ ^,. ?^ 306.... 



60 - 2i0iV4 • • - ^ • V 60 - 210 - 0.003 

Calculating the quadrature formula with A^ = 500, we obtain the approximate value 0.0062067..., 
whence (i'(4) > 0.0062067... - 2 - 0.003 > 0. D 

Lemma 15. We have d"{A) > 0. 

Remark 16. By numerical calculation, d"{A) ^ 0.0541341.... 

Proof. From ( IT9|) with t = 4,j = 2 and inserting the values of M^s given by ([7]) together 
with e < 3.7 we get ||i:^^^||oo < 7 - lO^^. 

Thus to bring the error below 6 = 0.027, we need to chose the step number N large enough 
to have 

^■^°'' <6 i.e. N>No:=(l^ l:!^[\._ ^255.... 



60 - 210 A^4 - V 60 - 210 - 0.027 

Calculating the quadrature formula with A^ = 500, we obtain the approximate value 0.05413417. 
whence d"{4) > 0.05413417... - 2 - 0.0027 > 0. D 

Lemma 17. We have (i"'(4) > 0. 

Remark 18. By numerical calculation, d'"{4) ^ 0.2255707.... 
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Proof. From (1191) with t = 4, j = 3 and calculating with the same values of M^ from ([7]) as 
above - together with £ < 3.7 in view of (1431) - we arrive at 

^iv ^ Qfif^QQ + 210L + 150^2 + 24L3) + QGC^CiQ + 54L + TSL^ + 2AL^) 

+ 4G^G'G"'{6L + 21L^ + 12L^) + G^G^^{?,L^ + AL^) + 3G^G"^(6L + 21L^ + 12L^), 

\H^^\ < 5.4- 10^1 

To bring this down below 5 = 0.112, we need to chose the step number A^ large enough to 
have 

'■'■">" <^ i.e. 7V>7Vo:=,'/I^1^^298. 



60-2i0iV4 • • - ^ • V 60 -210 -0.112 

Calculating the quadrature formula with A^ = 500, we obtain the approximate value 0.22557089. 
whence d'"(4) > 0.22557089... - 2 ■ 0.112 > 0. D 

So we arrive at the analysis of dX . Numerical tabulation of values give that dY is decreasing 
from (i^(4) f» —2,217868... to even more negative values as t increases from 4 to 5. So 
we now set forth proving that d^ < in [4,5]. To arrive at it, our approach will be 
a computation of some approximating polynomial p{t), which is, within a small and well 
controlled error, will be a Taylor polynomial of dy{t). However, as we intend to keep the 
step number A^ of the numerical integration under 500, we take the liberty of approximating 
d^ by different polynomials (using different Taylor expansions) on various subintervals of 
[4,5]. More precisely, we divide the interval [4,5] into 2 parts, and construct approximating 
Taylor polynomials around 4.25 and 4.75. 

So now setting to = 4.25 or to = 4.75, the Taylor approximation will have the form 

(45) ci^(t) = V A^(t-to)^ + i?„(rf^to,t), i?„(rf^to,t) = ^— ^(t_t^)"+i. 

Therefore instead of fl36l) we can use 



\l^ (^y + ^M ^ ll-^€,n+6,+ l|LM0,l/2] + ||-^g,«+6,-|UM0,l/2] .-(n+1) 

\Kn{a ,to,tj < — — 4 

(n + Ij! 



(46) < 



< 



i II W II _l_ i II ff II 

2 II ?i'^+6,+ ||oo T 2 II Cii+6,— ||oo 

(n + l)!22"+2 

KiaX|^_jg|<l/4 ||ifg^„+6,+ ||oo + maX|^_tp|<i/4 \^H^^n+%,- 



(n + 1)122^^+3 

So once again we need to maximize f lT2|) . that is functions of the type f^| logfl™, on [0,9] 
(or, more precisely, on the subinterval TZ{G) ~ [0.02776..., 9], where the values are actually 
attained by v := G{x)). So now similarly to fl39|) . we get from flT^ of Lemma |6]that for any 

m < 31 and |^ - to| < 1/4 

\\H^,m,±\\oo = max I (^) , 9« log" 9! = 9« log™ 9 < 9*°+^/^ log'" 9. 



In all, for n < 25 

.,7. Ilf, . .,, ./9^-^log"+^9 = 19,683 2"+6log'^+^3 if to = 4.25, 

^ ^ IC-'Sf#i/4 ll^«.-+6,±(Xj||oo S < 95lQgn+69 ^ 59, 049 2"+^ log"+6 3 if^^ ^ 4 75, 



In case to = 4.25 now we chose n = 7. Then for this case the Lagrange remainder term (H6l) 

SI 4 928 lop'"^^ 3 
of the Taylor formula fl^5l) can be estimated as \Rn(d^ ,t)\ < '—. < 0.21 =: 6s. 

V! — V I "A 71 - 2"(n + 1)! 



HARDY-LITTLEWOOD MAJORANT PROBLEM 17 

As before, the Taylor coefficients dj^^{tQ) cannot be obtained exactly, but only with some 
error, due to the necessity of some kind of numerical integration in the computation of the 
formula (l27j) . Hence we must set the partial errors ^o, • • • , ^7 with X]7=o ^i < ^ '■= 2.21, say, 
so that d^{t) < Pnit) + 5 for 



(48) Pn{t):='£^{t 

The analogous criteria to (1571) now has the form: 



4.25)^ 



(49) 



rf(j+5)(4.25)-rfj- 



(t - 4.25)^ 






<6j (j = 0,l,...,n). 



That the termwise error f H9|) would not exceed 6j will be guaranteed by Nj step quadrature 
approximation of the two integrals in f l27p defining d'^^~^^\4:.25) with prescribed error t]j each. 
Therefore, we set rij := 6j2'^^j\/2, and note that in order to have 



(50) 



iv, > iv; 



I H^"^ II 



60-2 



10, 



Vj 



I H^"^ II 

|-"4.25,j+5,±l|oo 

60 ■ 2i0j!22i"i5,- 



riv 



suffices by the integral formula fITT]) and Lemma O That is, we must estimate ||-f^4 25 0+5 ±11 00 
for j = 0, . . . , 7 and thus find appropriate values of N*. 

Lemma 19. For j = 0, ... ,7 we have the numerical estimates of Table for the values of 
ll-^4.2^5j+5,±lloo- Setting 6j as seen in the table for j = 0, . . . ,7, the approximate quadrature 
of order 500 := Nj > N* with the listed values of N* yield the approximate values dj 



as listed in Table admitting the error estimates fj49|) for j = 0, . . . , 7. Furthermore, 
||-R8(c^^;^)||oo < 0.21 =: 6s and thus with the approximate Taylor polynomial P^lt) defined in 
the approximation \d^{t) — P7{t)\ < 6 := 2.21 holds uniformly for t G [4,4.5]. 



Table 2. Estimates for values of 
suiting Nj, and dj with Nj := N := 



l-^4.2^5j+5,±lloo, chosen values of 6j and re- 
: 500 for j = 0, . . . , 7. 




1 

2 
3 

4 

5 
6 
1 



estimate for 
1.23 
5.32 
2.29 
9.80 
4.18 
1.77 
7.47 
3.14 



/7^^ I 

-"4.25,i+5.±l 

10^5 
10^5 
10^6 

10^9 



6, 


n; 


0.65 


499 


0.73 


494 


0.4 


492 


0.15 


486 


0.04 


486 


0.01 


466 


0.01 


302 


0.01 


188 



dj 
-11.99030682 
-64.72801527 
-273.5687453 
-1000.494741 
-3319.462864 
-10,266.25853 
-30,113.02268 
-84,761.00164 



Proof. We start with the numerical upper estimation of HI25 j ±ix) for A < x < 4.5. In the 



general formula fl20|) now we consider the case t = 4.25 and use the estimates ([7]) of Mi, M2, 
M3 and M4, together with i < 3.7 - c.f. fH5]) - to compute 

|-"4.25j,±(^) 



< 3.7-'|4.78 ■ lO^j^ + 3.72 ■ 10*^j^ + 1.09 ■ lO^^j^ + 1.44 ■ lO^^j + 7.29 ■ 10"|. 
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Finally, we collect the resulting numerical estimates of ||ii/^'^^|| in Table |2] and list the cor- 
responding values of A^* and dj, too, as given by the formulae (150|) and the numerical 
quadrature formula (TTTl) with step size h = 0.001, i.e. N = Nj = 500 steps. D 

Lemma 20. We have d^{t) < for all A < t < 4.5. 

Proof. We approximate d^{t) by the polynomial Prit) constructed in fHHl) as the approximate 
value of the order 7 Taylor polynomial of d^ around to '■= 4.25. As the error is at most 6, it 
suffices to show that p{t) := P7{t)+6 < in [4, 4.5]. Now Pj{4) = -2.2178666857... so P7(4) + 

6 <0. Moreover, p'{t) = P^{t) = ^J^^ . . \. (t-4.25y-^ andp'(4) = -20.41147631... < 0. 

U - 1)' 
From the explicit formula oip{t) we consecutively compute also p"{4) = —104.6546745... < 0, 

p'"(4) = -426.8260106... < 0, p(^)(4) = -1473.198415... < 0, p(5)(4) = -5386.784165... < 

and p(^)(4) = -8922.772271... < 0. Finally, we arrive at p^'^\t) = d7=-84,761. 00164... . We 

have already checked that ^'^■'^(4) < for j = ... 6, so in order to conclude p{t) < for 

4 < t < 4.5 it suffices to show p^'^\t) < in the given interval. However, p^"^^ is constant rfy, 

hence p^'^\t) < for all t e'K. It follows that also p{t) < for all t > 4. D 

In case of to = 4.75 we have for all ^ G [4.5, 5] 
(51) max I ('^') ,9«log"^9| <max|('-^)",9^2'"log"3} =9Mog"9 (Vm < 37). 

In all, \\H^^n+6,±{^)\\oo < 59,049 2"+'5log"+^3 for all 4.5 < ^ < 5 and 4 < n < 31. In view of 

944 784 loe'^"^^ 3 

im this yields \RJd^,t)\ < '— ^- < 9.1 =: 87 for n = 6. 

2"(n + 1)! 



(52) 



Next we set 5o, . . . ,Sq. Now the criteria ( 137|) is modified as 
d(^+^\A.75) -dj 



(t - 4.75)^ 



' ^ ^ -^<S, (j = 0,l,...,n) 



22^- J ! 

Since the numerical calculation gives that d^{4.5) ~ —39.96194643..., now we wish to chose 
the partial errors 5j so that J2]=o ^i < ^ '■= 39.9, say, so that d^ {t) < Pn{t) + 5 with 

" -1 
(53) P„(t):=^-f(t-4.75)^ 

j=o ■'' 

The goal is that the termwise error (l52l) would not exceed 5j, which will be guaranteed by 
Nj step quadrature approximation of the two integrals defining d^^^^\A.7b) with prescribed 
error rjj each. Therefore, we set rjj := 5j2'^^ j\/2., and note that in order to have flS2]) 



(54) jy. ^ ^v. ._ 4/ J^i75,j+5,±lloo _ 4/ ll^i.75,i+5,±lloo 



^ ^' ■ y 60 ■ 21%- y 60 ■ 2i0j!22i-i5j 



TlV 



suffices by the integral formula fITT]) and Lemma [51 That is, we must estimate ||-f^4.75j4.5,±iioo 
for j = 0, . . . , 7 and thus find appropriate values of N* . 

Lemma 21. For j = 0,...,6 we have the numerical estimates of Tahle\^for the values 
of ||-f^4.V5.j,±lloo- Setting 6j as can be seen in the table for j = 0, . . . , 6, the approximate 
quadrature of order 500 '■= Nj > N* with the listed values of Nj yield the approximate values 
dj as listed in Table\^ admitting the error estimates (152]) for j = 0,...,6. Furthermore, 
\\Rf{d^ ,t)\\oo < 9.1 =: 67 and thus with the approximate Taylor polynomial Peit) defined in 
(15 3 p the approximation \d^{t) — Peit)\ < 6 := 39.9 holds uniformly for t G [4.5, 5]. 
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-^4.75j+5,±ll«3' ^^^ values of 6j and the 



resulting N*, and the values of dj with N := Nj := 500 steps for j = 0, . . . , 6. 




1 

2 
3 

4 
5 
6 



estimate for 
4.98 


ttIV 

-"4.75j+5,± 

10^5 


2.13 


10^6 


9.07 


10^6 


3.85 


10^7 


1.63 


10^8 


6.83 


w 


2.86 


w 



6, 
8 


N* 
318 


9 


313 


7 


339 


3 


323 


1 


305 


1 


207 


1 


134 



dj 
-111.5230149 
-432.5730841 
-1509.259871 
-4867.920658 
-14,185.12009 
-42,842.09045 
-119,563.5221 



Proof. We start with the numerical upper estimation of Hl\^ , ^{x) for 4.5 < x < 5. In fl20l) 
now we insert t = 4.75, use again the estimates (jTj) of Mi — M4 and i < 3.7 and arrive at 



|-"4.75j,±(^) 



< 3.7-'|l.44 ■ lO^j^ + 1.23 ■ 10^ f + 3.93 ■ lO^^j^ + 5.7 • 10"j + 3.18 • lO^^j 



Finally, we collect the resulting numerical estimates of ||-ff^^|| in Tableland list the corre- 
sponding values of A^* and dj, too, as given by formulae (15^ and the numerical quadrature 
formula ([II]) with step size h = 0.001, i.e. N = Nj = 500 steps. D 

Lemma 22. We have d^{t) < for all 4.5 < t < 5. 

Proof. We approximate d^(t) by the polynomial Peit) constructed in (153|) as the approx- 
imate value of the order 6 Taylor polynomial of d^ around to •= 4.75. As the error 
is at most 6 = 39.9, it suffices to show that p(t) := Peit) + S < in [4.5,5]. Now 



Pg(4.5) = -39.9655627058... so Pg{4:.5) + 6 < 0. Moreover, p'{t) = P^{t) 



di 



l^j=i ^j _ 1) 



it- 



4.75)-' ^ and p'(4.5) = —174.8777051... < 0. From the explicit formula of p{t) we con- 
secutively compute also p"(4.5) = -662.2069802... < 0, p'"(4.5) = -2199.092624... < 0, 
p(4)(4.5) = -7810.957541. ..2 < and p^^\4:.5) = -12,951.20993... < 0. Finally, we arrive 
at p^^\t) = rf6=-119,563.5221... We have already checked that p(^)(4.5) < for j = . . . 5, 
so in order to conclude p{t) < for 4.5 < t < 5 it suffices to show p^^\t) < in the given 
interval. However, p^^' is constant, so p^^'{t) < for all t G M. It follows that also p{t) < 
for all t > 4.5. D 

5. Conclusion 

With the help of the sharper quadrature formula fITT]) further numerical analysis is possible 
for higher values of k. In principle we can divide the interval {k, k-\-l) to smaller and smaller 
intervals to get improved error estimations of Taylor expansions to compensate the larger 
and larger error bounds resulting from e.g. ([7]) and the increase of t. We have a strong feeling 
that this way we could work further to higher values of k. However, even that possibility 
does not mean that we would have a clear theoretical reason, a firm grasp of the underlying 
law, rooted in the nature of the question, for what the result should hold for all k. 
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